arXiv:1505.00034v2 [astro-ph.CO] 13 Sep 2015 
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We propose a heuristic unification of dark matter and dark energy in terms of a single “dark 
fluid” with a logotropic equation of state P = A\n(p/pp), where p is the rest-mass density, pp = 

5.16 x 10 99 gm -3 is the Planck density, and A is the logotropic temperature. The energy density 
e is the sum of a rest-mass energy term pc? mimicking dark matter and an internal energy term 
u(p) = —P(p) — A mimicking dark energy. The logotropic temperature is approximately given by 
A ~ pac 2 / \n(pp /pa) ~ pac 2 /[1231ii(10)], where pA = 6.72 x 10 ~ 24 gm -3 is the cosmological density 
and 123 is the famous number appearing in the ratio pp/pA ~ 10 123 between the Planck density 
and the cosmological density. More precisely, we obtain A = 2.13 x 10 ~ 9 gm -1 s -2 that we interpret 
as a fundamental constant. At the cosmological scale, this model fullfills the same observational 
constraints as the ACDM model (they will differ in about 25 Gyrs when the logotropic universe 
becomes phantom). However, the logotropic dark fluid has a nonzero speed of sound and a nonzero 
Jeans length which, at the beginning of the matter era, is about A j = 40.4 pc, in agreement with 
the minimum size of the dark matter halos observed in the universe. At the galactic scale, the 
logotropic pressure balances gravitational attraction and solves the cusp problem and the missing 
satellite problem. The logotropic equation of state generates a universal rotation curve that agrees 
with the empirical Burkert profile of dark matter halos up to the halo radius. In addition, it implies 
that all the dark matter halos have the same surface density Eo = porh = 141 Af 0 /pc 2 and that 
the mass of dwarf galaxies enclosed within a sphere of fixed radius r u = 300 pc has the same value 
A /300 = 1-93 x 10' M 0 , in remarkable agreement with the observations. It also implies the Tully- 
Fisher relation Mb/v\ = 44A/ 0 km _ 4 s 4 . We stress that there is no free parameter in our model 
(we predict the values of Eo, A /300 and A/f,/v 4 in terms of fundamental constants). We sketch 
a justification of the logotropic equation of state in relation to the Cardassian model (motivated 
by the existence of extra-dimensions) and in relation to Tsallis generalized thermodynamics. We 
also develop a scalar field theory based on a Gross-Pitaevskii equation with an inverted quadratic 
potential, or on a Klein-Gordon equation with a logarithmic potential. 

PACS numbers: 95.30.Sf, 95.35.+d, 95.36.+X, 98.62.Gq, 98.80.-k 


I. INTRODUCTION 

The nature of dark matter (DM) and dark energy (DE) 
is still unknown and remains one of the greatest myster¬ 
ies of modern cosmology. DM has been introduced in as¬ 
trophysics to account for the missing mass of the galax¬ 
ies inferred from the virial theorem [J and to explain 
their flat rotation curves [2]. DE has been introduced 
in cosmology to account for the present acceleration of 
the expansion of the universe [3]. In the standard cold 
dark matter (ACDM) model, DM is represented by a 
pressureless fluid and DE is ascribed to the cosmologi¬ 
cal constant A introduced by Einstein J]. The ACDM 
model works remarkably well at the cosmological scale 
but it encounters serious problems at the galactic scale. 
In particular, it predicts that DM halos should be cuspy 
[5] while observations reveal that they have a flat core [6]. 
On the other hand, the ACDM model predicts an over¬ 
abundance of small-scale structures (subhalos/satcllites), 
much more than what is observed around the Milky Way 
PFj. These problems are referred to as the “cusp prob¬ 
lem” and “missing satellite problem”. The expression 
“small-scale crisis of CDM” has been coined. 

There are also unexplained important observational re¬ 
sults. For example, it is an empirical fact that the sur¬ 
face density of galaxies has the same value Eo = po r h = 
141 ±f 2 A/q/pc 2 even if their sizes and masses vary by sev¬ 


eral orders of magnitude (up to 14 orders of magnitude 
in luminosity) 0- On th e other hand, it is known that 
the asymptotic circular velocity of the galaxies is related 
to their baryonic mass by the Tully-Fisher (TF) relation 
Mb/ V h = 47 ± 6Af0km _4 s 4 [3] HU] , Finally, Strigari et 
al. nn have shown that all dwarf spheroidal galaxies 
(dSphs) of the Milky Way have the same total DM mass 
contained within a radius of r u = 300 pc. From the obser¬ 
vations, they obtained log(A/ 30 o/Af 0 ) = T.O^q^. To our 
knowledge, there is no theoretical explanation of these 
observational results. 

The small scale problems of the ACDM model are re¬ 
lated to the assumption that DM is pressureless. This 
assumption is valid if DM is made of weakly interacting 
massive particles (WIMPs) with a mass in the GeV-TeV 
range. These particles freeze out from thermal equilib¬ 
rium in the early universe and, as a consequence of this 
decoupling, cool off rapidly as the universe expands. In 
order to solve the small-scale crisis of CDM, some authors 
have developed alternative models of DM. For example, 
it has been proposed that DM halos are made of fermions 
(such as sterile neutrinos) with a mass in the keV range 
mug, or bosons (such as axions) in the form of Bose- 
Einstein condensates (BECs) with a mass ranging from 
10~ 2 eV to 10 -20 eV depending whether the bosons inter¬ 
act or not [14]. In these models, the quantum pressure 
prevents gravitational collapse and leads to cores instead 
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of cusps. These models sometimes provide a good fit of 
the rotation curves of galaxies but they do not explain 
the universality (and the values) of E 0 , M^/v^, and M 300 . 

On the other hand, at the cosmological scale, despite 
its success at explaining many observations, the ACDM 
model has to face two theoretical problems. The first one 
is the cosmic coincidence problem, namely why the ratio 
of DE and DM is of order unity today if they are two 
different entities m- The second one is the cosmologi¬ 
cal constant problem [16j. The cosmological constant A 
is equivalent to a constant energy density ca = Pac 2 = 
Ac 2 /87tG associated with an equation of state P = — e 
involving a negative pressure. Some authors P3 have 
proposed to interpret the cosmological constant in terms 
of the vacuum energy. Cosmological observations lead to 
the value pA = A/ 87 rG = 6.72 x 10 -24 g nr -3 of the cosmo¬ 
logical density (DE). However, particle physics and quan¬ 
tum field theory predict that the vacuum energy should 
be of the order of the Planck density pp = c 5 /HG 2 = 
5.16 x 10"gm' 3 . The ratio between the Planck density 
pp and the cosmological density pa is 

^ ~ 10 123 , ( 1 ) 

PA 

so these quantities differ by 123 orders of magnitude! 
This is the origin of the cosmological constant problem. 
To circumvent this problem, some authors have proposed 
to abandon the cosmological constant A and to explain 
the acceleration of the universe in terms of a dark en¬ 
ergy with a time-varying density associated with a scalar 
field called “quintessence” [TS]. As an alternative to 
quintessence, Kamenshchik et al. ns have proposed a 
heuristic unification of DM and DE in terms of an ex¬ 
otic fluid with an equation of state P = —A/e called the 
Chaplygin gas. This equation of state provides a model 
of universe that behaves as a pressureless fluid (DM) at 
early times, and as a fluid with a constant energy density 
(DE) at late times, yielding an exponential acceleration 
similar to the effect of the cosmological constant. How¬ 
ever, in the intermediate regime of interest, this model 
does not give a good agreement with the observations 
12Uj so that various generalizations of the Chaplygin gas 
model have been considered. In this Letter, we propose 
a new model based on a logotropic equation of state |21] 
that seems to give a solution to all the problems men¬ 
tioned above and, most importantly, that predicts the 
correct values of E 0 , M^/v^, and M 300 with remarkable 
accuracy, and without free parameter. 


II. LOGOTROPIC COSMOLOGY 


A. The logotropic dark fluid 


The Friedmann equations for a flat universe without 
cosmological constant are [22]: 


de 

dt 


+ 3-(e + P) =0, 
a 


H 2 


/a \ 2 87 tG 
\ay 3c 2 £ ’ 


( 2 ) 


where e(t) is the energy density, P(t) is the pressure, a(t) 
is the scale factor, and H = a/a is the Hubble parameter. 

For a relativistic fluid at T = 0, or for an adiabatic 
evolution (which is the case for a perfect fluid), the first 
law of thermodynamics reduces to j22| : 


de = 


P + e 
P 


dp , 


( 3 ) 


where p is the rest-mass density. Combined with the 
equation of continuity we get 


dp 

dt 


o a 

3 -p 
a 


0 => p 


Po 


( 4 ) 


where p 3 is the present value of the rest-mass density, and 
the present value of the scale factor is taken to be a 3 = 1. 
This equation, which expresses the conservation of the 
rest-mass, is valid for an arbitrary equation of state. 

For an equation of state specified under the form P = 
P(P), Eq. © can be integrated to obtain the relation 
between the energy density e and the rest-mass density. 
We obtain 


e = pc~ + p 


dp' = pc 2 + u(p), (5) 


where the constant of integration is set equal to zero. We 
note that u(p) can be interpreted as an internal energy 
density. Therefore, the energy density e is the sum of 
the rest-mass energy pc 2 and the internal energy u{p). 
The rest-mass energy is positive while the internal energy 
can be positive or negative. Of course, the total energy 
e = pc 2 + u(p) is always positive. 

We assume that the universe is filled with a single dark 
fluid described by a logotropic equation of state 


p=m "{£)- (6) 

where A is the logotropic temperature (determined be¬ 
low) and pp = 5.16 x 10" gm -3 is the Planck density. 
It will be called the Logotropic Dark Fluid (LDF). Using 
Eqs. © and ©, the relation between the energy density 
and the rest-mass density is 


e = pc 2 — Ain ^ — A = pc 2 +u(p). (7) 

The energy density is the sum of two terms: a rest-mass 
energy term pc 2 oc a~ 3 that mimics DM and an internal 
energy term u(p) = —P(p) — A that mimics DE. This de¬ 
composition leads to a natural, and physical, unification 
of DM and DE and elucidates their mysterious nature. 
We note that the pressure is related to the internal energy 
by P = —u — A. Combining Eqs. © and ©, we obtain 
e = ppc 2 e p / A — P — A which determines, by inversion, the 
equation of state P(e). From Eqs. ©, © and ©, we get 
P = A\n(po/ppa 3 ) and e = p 3 c 2 /a — Aln(po/ P pa 3 ) — A. 
We note that the internal energy u = — A\n(p/pp) — A 
is positive for p < pp/e and negative for p > pp/e. 
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In the early universe (a —► 0, p —► +oo), the rest-mass 
energy (DM) dominates, so that 

P~Alnf— ( 8 ) 

a 3 \Ppc z J 

For small values of the scale factor, we recover the results 
of the CDM model (P = 0) since e oc a -3 . In the late 
universe (a —> + 00 , p —>• 0), the internal energy (DE) 
dominates, and we have 


density pp and the cosmological density p A (itself ob¬ 
tained from the Hubble constant Hq and the fraction of 
DM fl m0 )- From now on, we shall regard A as a fun¬ 
damental constant that supersedes the cosmological con¬ 
stant. We note that it depends on all the fundamental 
constants of physics H, G, c, and A [see Eqs. (fTTb and 

© • 

After simple manipulations, the rest-mass density, the 
pressure and the energy density of the LDF can be ex¬ 
pressed in terms of B as 


e ~ —Ain ^ —^ ~ 3Alna, P ~ —e. (9) 

We note that the equation of state P(e) behaves asymp¬ 
totically as P ~ — e, similarly to the usual equation of 
state of DE. It is interesting to recover the equation of 
state P = — e from the logotropic model ©• This was 
not obvious a priori. 


B. The logotropic temperature 


pc 2 _ 0 

e 0 a 6 


P 

£A 


-B-l + Bln 


PC 


\ CqH 


m, 0 


, (13) 


— =-B- l-3Hlna, (14) 

£A 


e 


eo 


“-1“ (1 ~ ^ra,o) 

eo 


1 + B In 




;,0 6 0 


pc z 


(15) 


Since, in our model, the rest-mass energy of the dark 
fluid mimics DM, we identify po with the present density 
of DM. We thus set p 0 = n m , 0 e 0 /c 2 = 2.54 x 10 -24 gm -3 , 
where e 0 /c 2 = 3Hf/8nG = 9.26 x 10~ 24 gm- 3 is the 
present energy density of the universe (we have taken 
H 0 = 70.2 km s -1 Mpc” 1 = 2.275 10" 18 s" 1 ) and fi m , 0 = 
0.274 is the present fraction of DM (we also include 
baryonic matter). As a result, the present internal en¬ 
ergy of the dark fluid uq/c 2 = eo/c 2 — po is identified 
with the present density of DE p A = (1 — fl mi0 )e 0 /c 2 = 
6.72 x 10~ 24 g m -3 where Ha,o = 1 — fi m ,o = 0.726 is the 
present fraction of DE. 

Applying Eq. 0 at a = 1, we obtain the identity 

PP = O 0 el+1/ B (10 ) 

Pa 1 


where we have defined the dimensionless logotropic tem¬ 
perature B through the relation A = Bp^c 2 . This iden¬ 
tity is strikingly similar to Eq. ([l]) which appears in rela¬ 
tion to the cosmological constant problem. In the present 
context, the identity (101 determines the logotropic 
temperature B. Qualitatively, B ~ l/ln(pp/pA) — 
l/[1231n(10)[. This gives a new interpretation to the 
famous number 123 ~ log(pp/pA) as being the inverse 
logotropic temperature. More precisely, we obtain 


B = 


1 


In (TAWiap") _ 1 
\ “m,0 PA J 


= 3.53 x 10 


-3 


(ii) 


and 


A = Bp A c 2 = 2.13 x 10' 


1 gm 1 s 2 . 


( 12 ) 


As a result, there is no free parameter in our model. The 
logotropic temperature is determined from the Planck 


e 

eo 


a 


m, 0 


+ (l-fim,o)(l + 3Blno), (16) 


- = D m , 0 e (B+1)/B e p / BeA -(l-H m ,o) (— + B). (17) 

e 0 \ e A / 

The ACDM model is recovered for B = 0, i.e., e/eo = 
pc 2 /e 0 + (1 - e/eo = fl m ,o/a 3 + 1 - and 

P = —ca- The ACDM model is equivalent to a con¬ 
stant negative pressure P = — e A H2) and to the relation 
e = pc 2 +e A between the energy density and the rest-mass 
density. According to Eq. (10), the condition B = 0 in 


the logotropic model corresponds to pp = + 00 , hence 
h = 0. Therefore, the fact that B is small but nonzero 
as vindicated by the observations (see below) shows that 
quantum mechanics (h ^ 0) plays a role in the late uni¬ 
verse in relation to DE. 


C. Evolution of the logotropic universe 


The relation between the energy density and the rest- 
mass density [see Eq. is plotted in Fig. [lj The 

evolution of the energy density with the scale factor [see 
Eq. (16)] is plotted in Fig. [ 2 ] The universe starts at 
a = 0 with an infinite rest-mass density (p —¥ + 00 ) and 
an infinite energy density (e —> + 00 ). 1 The rest-mass 
density decreases as a increases. The energy density first 


1 Of course, our model that attemps to unify DM and DE is only 
valid at sufficiently late times, typically for a > a; = 10 4 , after 
the inflation and the radiation eras. Therefore, the limit a —> 0 
is here formal. 
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decreases as a increases (i.e. p decreases), reaches a min¬ 
imum cm = — A\n(A/ppc 2 ) at a« = (poc 2 /A) 1 / 3 (i.e. 
Pm = A/c 2 ), then increases as a increases (i.e. p de¬ 
creases) further, and tends to e +oo as a —> +oo (i.e. 
p —> 0). The branch a < cim (he. p > pm) corresponds 
to a normal behavior in which the energy density de¬ 
creases as the scale factor increases. The branch a > aM 
(i.e. p < pm) corresponds to a phantom behavior [24j in 
which the energy density increases as the scale factor in¬ 
creases. We note that A is equal to the rest-mass energy 
at the point where the universe becomes phantom. 



FIG. 1: Relation between the energy density e and the rest- 
mass density p in the logotropic model. It is compared with 
the relation e = pc 2 + ea corresponding to the ACDM model. 
The energy density presents a minimum (e/eo )m = 0.7405 at 
pM<?/t o = 2.56 x 10 -3 separating the normal universe and 
the phantom universe. 



FIG. 2: Evolution of the energy density as a function of the 
scale factor in the logotropic model. It is compared with 
the ACDM model. The energy density presents a minimum 
(e/eo )m = 0.7405 at aju = 4.75. 

The evolution of the pressure with the scale factor [see 
Eq. @] is plotted in Fig. [3j The pressure decreases as 



FIG. 3: Evolution of the pressure as a function of the scale fac¬ 
tor in the logotropic model. It is compared with the ACDM 
model where P = — ea- The pressure becomes negative at 
a w = 7.00 x 10 -42 The point separating the normal uni¬ 
verse from the phantom universe is located at hm = 4.75 
and Pm/ t a = —1 .02. 



FIG. 4: Equation of state P(e) of the logotropic model. It is 
compared with the equation of state P = — ea of the ACDM 
model. 


a increases (i.e. p decreases). It starts from P -T +oo 
at a = 0 (i.e. p —> +oo, e —> +oo), vanishes at a w = 
(po/pp) 1 / 3 , achieves the value Pm = at aM (he. 
Pm, cm) and tends to P — > — oo as a — > +oo (i.e. p —> 0, 
e —> Too). The equation of state P(e) [see Eq. @] is 
defined for e > 6 m and has two branches corresponding 
to a normal universe {P > Pm) and a phantom universe 
(P< Pm), as shown in Fig. [4] Therefore, the equation 
of state P(e) is multi-valued. 

The speed of sound c s , defined by c 2 s = P'(e)c 2 , is 
given by c 2 /c 2 = 1 /{pc 2 /A - 1) = l/[(a M /a) 3 - 1]- It 
is real for a < aM (i.e. when the universe is normal) 
and imaginary for a > «« (i.e. when the universe is 
phantom). The relation between the speed of sound and 
the scale factor is plotted in Fig. [5j 

Solving the Friedmann equation Q with the energy 
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FIG. 5: Evolution of the speed of sound with the scale factor 
in the logotropic model. For the ACDM model, c a = 0. The 
speed of sound is equal to the speed of light (c s = c) at 
as = aM/2 1 / 3 = 3.77. At the present time (a = ao = 1), we 
have ( c s /c ) 2 = 1 /(a 3 M — 1) = 9.42 x 10~ 3 . 



FIG. 6: Temporal evolution of the scale factor in the lo¬ 
gotropic model as compared to the ACDM model. The 
age of the universe is to = 13.8 Gyrs. In the logotropic 
model, the universe accelerates at t c = 7.19 Gyrs, the speed 
of sound exceeds the speed of light at t s = 34.5 Gyrs, and 
the universe becomes phantom at tM = 38.3 Gyrs (at present 
wo = Po/to = —0.729). 


density (16), we obtain the evolution of the scale fac¬ 
tor represented in Fig. [6] For t 0, we have 
a ~ (3^/Q m:0 H 0 t/2) 2 / 3 and e/e 0 ~ 4/(9 Hgt 2 ) as in the 
Einstein-de Sitter (EdS) model. For t —> +oo, we have 
a oc e 3m-tim,o)H 2 a P/4, and e / €o „ [3£ (1 _ Q mfi )H 0 t/ 2] 2 . 
This asymptotic solution, which is valid in the regime 
where the universe is phantom, has a super-de Sitter be¬ 
havior. There is no future finite time singularity (no Big 
Rip). The scale factor, the energy density and the pres¬ 
sure become infinite in infinite time (Little Rip). 


D. Cosmological implications 


From the observational viewpoint, there is no visible 
difference between the logotropic model and the ACDM 
model at large scales. Differences will appear in about 
25 Gyrs, when the universe becomes phantom (this as¬ 
pect will be developed in a future work). However, this 
moment is very remote in the future, and for the time be¬ 
ing, the logotropic model and the ACDM model behave 
similarly (see Figs. [2] and [6]) . This is satisfactory since 
the ACDM model works well at the cosmological scale. 
On a theoretical point of view, the logotropic model has 
several advantages with respect to the ACDM model. In 
our model, DM and DE are the manifestation of a sin¬ 
gle dark fluid described by a unique equation of state. 
Therefore, there is no cosmic coincidence problem. On 
the other hand, the cosmological constant problem of Eq. 
(jT|) is translated into an equation (10) that determines the 
logotropic temperature A ~ ca/[123 ln(10)]. 

An important difference between the ACDM model 
and the logotropic model concerns the speed of sound 
c s defined by c 2 /c 2 = P'(e). In the ACDM model, since 
P = 0 (actually, P = — ca), the speed of sound c s = 0. 
As a result, the Jeans length is zero (A j = 0), imply¬ 
ing that the homogeneous background is unstable at all 
scales so that halos of any size should be observed in 
principle. However, this is not the case. There does not 
seem to be halos with a size smaller than I? m i n ~ 10 pc. 
Contrary to the ACDM model, the logotropic model has 
a nonzero speed of sound, hence a nonzero Jeans length. 
We can obtain an estimate of the Jeans length A j at 
the beginning of the matter era where perturbations 
start to grow. We assume that the matter era starts 
at cti = 10 4 , corresponding to the epoch of matter- 
radiation equality. In this era, we can make the ap¬ 
proximation e = pc 2 , so the Jeans wavenumber is given 
by k 2 = Ai:Gpa 2 /c 2 [22], where c 2 = P'(p) = A/p. 
At a,i = 10 —4 , we find pi = 2.54 x 10 _12 g/m 3 and 
(c 2 /c 2 )i = 9.33 x 10~ 15 . This leads to a Jeans length 
A j = 2-K/kj = 1.25 x 10 18 m = 40.4pc which is of the 
order of magnitude of the smallest known dark matter ha¬ 
los such as Willman I (ry, = 33pc) (see Table 2 of f25|). 
We predict that there should not exist halos of smaller 
size since the perturbations are stable for A < A j. This 
is in agreement with the observations, unlike the ACDM 
model. Therefore, a small but finite value of B 1 yielding 
a nonzero speed of sound and a nonzero Jeans length, is 
able to account for the minimum observed size of dark 
matter halos in the universe. It also puts a cut-off in the 
density power spectrum of the perturbations and sharply 
suppresses small-scale linear power. This may be a way 
to solve the missing satellite problem. 


III. LOGOTROPIC DARK MATTER HALOS 

The interest of the logotropic model becomes evi¬ 
dent when it is applied to DM halos. We assume that 
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DM halos are described by the logotropic equation of 
state ^ with the logotropic temperature A = 2.13 x 
10~ 9 gm _1 s -2 determined previously, viewed as a funda¬ 
mental constant. At the galactic scale, we can use New¬ 
tonian gravity. Combining the condition of hydrostatic 
equilibrium VP + pV$ = 0 with the Poisson equation 
A<f> = 4nGp , assuming spherical symmetry, and intro¬ 
ducing the notations 6 = po/p and £ = r/ro, where po is 
the central density and 


r o = 


A 

4 


1/2 


(18) 


is the logotropic core radius, we obtain the differential 
equation 


1 d / 2 d0\ 1 

di) = ~e 


(19) 


with 0(0) = 1 and 0'(O) = 0. It can be viewed as a 
Lane-Emden equation of index n = — 1 [26j . This equa¬ 
tion has a simple analytical solution p s = (A/87rG) 1 / 2 r _1 
called the singular logotropic sphere because it diverges 
at the origin [21]. 2 The regular solutions must be com¬ 
puted numerically. They have a flat core and behave as 
p ~ (A/87rG) 1 / 2 r _1 for r —» + 00 . Since the logotropic 
spheres are homologous, they generate a universal DM 
profile. Indeed, if we rescale the density by the central 
density po and the radial distance by the core radius r 0 , 
we get an invariant density profile l/0(£). We note that 
the total mass of a logotropic sphere is infinite because 
of the slow decay of the density. This means that the lo¬ 
gotropic distribution cannot describe the whole cluster. 
It is valid only in the core. At larger distances, we must 
take into account complex physical processes such as tidal 
effects and incomplete relaxation that steepen the den¬ 
sity profile (see, e.g., [13]). It is usually found that the 
density profiles of DM halos decrease at large distances 
as r -3 0E]- Since we do not take these complicated 
processes into account, our logotropic model is only valid 
up to a few values of the core radius Vq. However, this 
is sufficient to determine the physical characteristics of 
DM halos. 

Using the Lane-Emden equation (19) the mass pro¬ 
file M(r) = fg p(r') Airr' 2 dr' is given by M[r) = 
47rpo r o£ 2 f? , (£)- The circular velocity defined by u 2 (r) = 
GM(r)/r can be expressed as u 2 (r) = AnGpor^O' {£,). 
We define the halo radius rh as the radius at which 
p/po = 1/4. The dimensionless halo radius is the so¬ 
lution of the equation 9(t; h ) = 4. We numerically find 
£ h = 5.8458 and 0’{£ h ) = 0.69343. Then, r h = £ h r 0 . The 
normalized halo mass at the halo radius is given by 


M h A 0'(6O 

- 3 = 4n —t - 

Po^ & 


1.49. 


( 20 ) 


2 We note, parenthetically, that this singular solution oc r 1 is 
similar to NFW cusps [5]. 


This value is relatively close to the value Mh/p$r\ = 1.60 
[T2l [13] obtained from the empirical Burkert profile [6] 
that provides a good fit of DM halos. On the other 
hand, the universal rotation curve predicted by the lo¬ 
gotropic model is very close to the Burkert profile up to 
the halo radius, i.e. for r < ru (see Fig. [?|. Very re¬ 
cently, Burkert [27] observed that the density profile of 
real DM halos behaves approximately as r _1 close to the 
halo radius. Interestingly, we note that the exponent —1 
precisely corresponds to the characteristic exponent of 
the logotropes. 



FIG. 7: Circular velocity profile of a logotropic sphere (solid 
line) in the region r < r^ where the model is valid. It is 
compared with the Burkert profile (dashed line). 


In addition to these already encouraging results, the 
logotropic equation of state has a very interesting prop¬ 
erty. According to Eq. (18), the surface density of the 


logotropic sphere is given by 

E » s «”■» = (sb) 


1/2 


61 • 


( 21 ) 


Since the logotropic temperature A is the same for all 
the halos (a consequence of our approach where we view 
A as a fundamental constant), this implies that the sur¬ 
face density of the DM halos should be the same. This 
is precisely what is observed [5 . Using the value of 
the logotropic temperature given by Eq. (12), we get 


Eq 11 = 141 Mq /pc 2 which coincides with the best-fit value 

-'obs 


Eq = 141-52 Ms/pc 2 °f the surface density of DM ha- 
los 0. This agreement is remarkable since there is no 
free parameter in our model. Furthermore, it is non triv¬ 
ial since the constant A depends, through Eqs. © and 
on the Planck density pp and on the cosmological 
density p\. This suggests that there is something deep 
behind these relations. 

There are interesting consequences of this result. Ac¬ 
cording to Eq. (20), the mass of the halos calculated 
at the halo radius rh is given by Mh = 1.49E 0 r|. On 
the other hand, the circular velocity at the halo ra¬ 
dius is u 2 = GMh/rh = 1.49E 0 Gr/j. Since the surface 
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density of the dark matter halos is constant, we ob¬ 
tain M/j/Mq = 210(rh,/pc) 2 oc r 2 and (wh/kms -1 ) 2 = 
0.905(r/j/pc) a r h . These scalings are consistent with the 
observations |12j . Furthermore, introducing the baryon 
mass Mb = fbMh where fb ~ 0.17 is the cosmic baryon 
fraction [10], we get 


Mb __ fb _ _ fb _ ,„o'| 

vi ~ 1.49E 0 G 2 ~ e'(£ h )(Ai tAG 3 ) 1 / 2 ' [ " ’ 

Therefore, v\ oc M/, which is the TF relation. More 
precisely, we predict (Mb/v^) th = 44MQkm _4 s 4 which 
is close to the observed value (Mb/vfy ohs = 47 ± 
6Mgkm~ 4 s 4 * * [TDj (we obtain a perfect agreement by tak¬ 
ing/ fe = 0.18). 3 * 

The logotropic equation of state also explains the 
observation of Strigari et al. El that all dSphs of 
the Milky Way have the same total dark matter mass 
M 300 contained within a radius r u = 300 pc, namely 
log(M^/M©) = 7.0+°; 3 . Using M 300 = ^p 0 roC e '(^u), 
= r u /r 0 and r 0 = r h /£h, we obtain M 300 = 
4nporh(rl/€h)6' (£hr u /rh). The logotropic distribution 
has the asymptotic behavior 0(£) ~ £/\/2 for £ —► +00 
m- For the dSphs considered in HH> f,hr u /r h » 1 (see 
Table 2 of [25]) so O'^hfu/^h) can be replaced by its 
asymptotic value 1/V2. This yields 


M 300 — 


4ttUo r 2 u 



(23) 


which is a constant, in agreement with the claim of Stri¬ 
gari et al. nn. We note that the constancy of M 300 is 
due to the universality of A. Furthermore, the numeri¬ 
cal application gives M 3 J} 0 = 1.93 x 10' M 0 , leading to 
log(M 3 Q 0 /M Q ) = 7.28 in very good agreement with the 
observational value. 

In conclusion, the logotropic equation of state can si¬ 
multaneously account for cosmological constraints (with 
the same level of precision as the ACDM model) and ex¬ 
plain properties of DM halos (their minimum size R m i n , 
their surface density E 0 , their mass M 30 0 , and the TF ra¬ 
tio Mb/v ^) that were not explained so far. This may be 
a hint that DM and DE are the manifestations of a single 
dark fluid. The best illustration of this “unification” is 
that we have obtained the value of A from cosmological 
constraints [see Eqs. © and ©] , and that this value 
accounts for the universality of the surface density Eg and 


3 The TF relation is sometimes justified by the MOND theory m 

which predicts a relation of the form = Ga^M^ between the 

asymptotic circular velocity and the baryon mass, where ao is a 

critical acceleration. Our results imply a|j = 1.72 x 10“ 10 m/s 2 

which is close to the value ag bs = (1.3 ± 0.3) x 10 _1 ° m/s 2 ob¬ 

tained from the observations Using the results of the Ap¬ 

pendix, we have aQ h ~ Hqc/A which explains why ao is of the 

order Hoc. We emphasize, however, that we do not use the 
MOND theory in our approach. 


mass M 30 o of DM halos [see Eqs. plj ) and (23)], as well 
as for the TF relation [see Eq. p2|]~~Assuming that this 
agreement is not a coincidence (the perfect agreement be¬ 
tween the predicted values of E 0 , M 3 oo, Mb/v £ and the 
observations is a strong support to our approach), the 
next step is to justify the logotropic equation of state. 
We sketch below several possible justifications based on 
extra-dimensions (Cardassian model), generalized ther¬ 
modynamics, and field theory. 


IV. POSSIBLE JUSTIFICATIONS OF THE 
LOGOTROPIC EQUATION OF STATE 

A. Cardassian model 


Freese and Lewis [25], in their so-called Cardassian 
model, have proposed to explain the accelerated expan¬ 
sion of the universe in terms of a modified Friedmann 
equation of the form 


H = 


8t tG 


'( P ), 


(24) 


where p = po a~ 3 is the rest-mass density and v(p) is a 
“new” term which characterizes the model. In the early 
universe, the term v{p) is negligible and one recovers the 
usual Friedmann equation of pressureless matter lead¬ 
ing to a decelerating universe with a oc t 2 / 3 (Einstein- 
de Sitter solution). In the late universe, the term zz(p) 
dominates and causes the universe to accelerate. Freese 
and Lewis [25! justify the modified Friedmann equation 


(24) in relation to the existence of extra-dimensions. The 


usual Friedmann equation is modified as a consequence 
of embedding our universe as a three-dimensional surface 
(3-brane) in higher dimensions. Our approach provides 
another, simpler, justification of this equation from the 
ordinary four dimensional Einstein equations. Starting 
from the usual Friedmann equation © and considering 
a dark fluid at T = 0, or an adiabatic fluid, with the 
energy density given by e = pc 2 + u(p) [see Eq. ©], we 
obtain Eq. (24) with u{p) = ( 8nG/3c 2 )u(p ). Therefore, 


the “new” term in the “modified” Friedmann equation 


(24) can be interpreted as the internal energy u(p) of 


the dark fluid while the “ordinary” term 8nGp/3 cor¬ 
responds to its rest-mass energy density pc 2 . There¬ 
fore, our approach provides a new justification of the 
Cardassian model. Inversely, the original justification 
of the Cardassian model, namely that the term v(p) 
arising in the modified Friedmann equation (24) may 


result from the existence of extra-dimensions, could be 
a way to justify the logotropic model corresponding to 
v(p) = -~{8 'kGA/8c 2 ) [In (p/pp) + 1]. In this respect, we 
note that the logotropic model asymptotically yields an 
equation of state of the form P ~ — e [see Eq. ([ 9 ])]. 
Using the virial theorem, one can easily show that this 
equation of state arises from a long-range confining force 
Fij = —3Uor 2 j that could be a fifth force or an effective 
description of higher dimensional physics [29]. 













B. Generalized thermodynamics 

The logotropic equation of state was introduced phe¬ 
nomenologically in astrophysics by McLaughlin and Pu- 
dritz p?fT j to describe the internal structure and the av¬ 
erage properties of molecular clouds and clumps. It was 
also studied by Chavanis and Sire [2T| in the context 
of Tsallis generalized thermodynamics m where it was 
shown to correspond to a polytropic equation of state of 
the form P = Kp 1 with 7 —> 0 and K —► 00 in such a way 
that A = 7 K is finite. It is associated with a generalized 
entropy of the form 


S L = / hi p dr , 


(25) 


which is called the Log-entropy m- The free energy can 
be written as iy = E — ASl, where E = ( 1 / 2 ) f p<f> dr is 
the gravitational energy. A critical point of iy at fixed 
mass M = f pdr 1 determined by the Euler-Lagrange 
equation SFl — p6M = 0, where p is a Lagrange mul¬ 
tiplier (chemical potential), leads to the Lorentzian-type 
distribution p(r) = l/[a + <b(r )/A\, where a = —p/A. 
We can check that this equation is equivalent to the 
equation of hydrostatic equilibrium with the logotropic 
equation of state (| 6 |. When combined with the Pois¬ 


son equation, we recover the Lane-Emden equation (19). 


These considerations show that A can be interpreted as 
a generalized temperature. This is why we call it the 
logotropic temperature. As a result, the universality of 
A (which explains the constant values of £ 0 , A /300 and 
Mb/vfr) may be interpreted by saying that the universe 
is “isothermal”, except that isothermality does not refer 
to a linear equation of state but to a logotropic equation 
of state in a generalized thermodynamical framework. If 
our model is correct, it would be a nice confirmation of 
the interest of generalized thermodynamics [31] in physics 
and astrophysics. 


C. Scalar field theory 

The logotropic equation of state can also be justified 
from a scalar field theory. If we view the dark fluid as a 
scalar field representing BECs, its evolution is described, 
in the nonrelativistic regime appropriate to DM halos, 
by the Gross-Pitaevskii equation [32]. The GP equation 
associated with the logotropic equation of state takes the 
form 


d'lb h~ 1 

ih m = - + - Am w'*’• (26) 

where ip(r,t) is the wave function. It is coupled to the 
Poisson equation A4> = 4ttG\iP\ 2 . Eq. (26) can be 
viewed as a GP equation with an inverted quadratic 
potential, i.e. with the exponent —2 instead of +2 in 
the usual GP equation [33]. To our knowledge, this 
equation has not been introduced before. Using the 


Madelung [34] transformation, this equation can be writ¬ 
ten in the form of fluid equations, called the quantum 
Euler equations [32], incorporating a quantum poten¬ 
tial Q = —(H 2 /2m)A^/~p/^/~p due to the kinetic term 
— (h 2 /2m)Aip (Heisenberg) and an isotropic pressure P 
due to the interaction term Am\ip\~ 2 ip . For the in¬ 
verted quadratic potential, the pressure P is given by 
the logotropic equation of state ©>• In the Thornas- 
Ferrni approximation, one can neglect the quantum po¬ 
tential. In that case, we recover the classical Euler- 
Poisson equations. A steady state of these equations 
satisfies the condition of hydrostatic equilibrium. Fur¬ 
thermore, one can show that the free energy Fl and the 
mass M are conserved by the Euler-Poisson equations. 
As a result, a minimum of free energy iy at fixed mass 
M is a stable steady state of the Euler-Poisson equa¬ 
tions [32]. This makes a correspondance with general¬ 
ized thermodynamics and shows the self-consistency of 
our approach. We also note that the GPP equations can 
be obtained as the nonrelativistic limit c —> +00 of the 
Klein-Gordon-Einstein equations [35]. The KG equation 
corresponding to the logotrope has a logarithmic poten¬ 
tial V(|0|) = — 2 Aln \(j>\ and writes 


1 d 2 <j) 
c 2 dt 2 


A (j> + 


m 2 c 2 

K 2 




0. (27) 


V. CONCLUSION 

We have proposed a heuristic unification of DM and 
DE in terms of a single dark fluid with a logotropic equa¬ 
tion of state (LDF). According to our model, what we 
usually call DM corresponds to the rest-mass density of 
the dark fluid and what we usually call DE corresponds 
to the internal energy of the dark fluid. 

At the cosmological scale, our model satisfies the same 
observational constraints as the ACDM model but avoids 
the cosmic coincidence problem (since DM and DE are 
the manifestation of a single dark fluid) and the cos¬ 
mological constant problem (since there is no cosmolog¬ 
ical constant in our approach). It also has a nonzero 
speed of sound and a nonzero Jeans length (contrary to 
the ACDM model) which can explain the minimum size 
-Rmin ~ 10 pc of DM halos. 

At the galactic scale, the logotropic pressure balances 
gravitational attraction and solves the cusp problem and 
the missing satellite problem of the CDM model. On the 
other hand, the logotropic model generates a universal 
rotation curve that provides a good agreement with the 
Burkert profile up to the halo radius. Furthermore, it 
implies that the surface density of DM halos and the 
mass of dwarf halos are the same for all the halos, in 
agreement with the observations. It also implies the TF 
relation. 

The most striking property of the logotropic model 
is the following. Using cosmological observations, we 
can obtain the value of the logotropic temperature A = 
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2.13 x 10 -9 gm -1 s -2 [see Eq. (12)]. It may be viewed 
as a fundamental constant since it actually depends on 
all the fundamental constants of physics H, G, c, and 
A. Then, applying the logotropic model to DM halos, 
and using thi s value of A, we ca n ob tain the value of So 
[see Eq. (21)], M b /v b [see Eq. (22)] and M 30 o [see Eq. 
(23)] which are in perfect agreement with the observa¬ 
tions. Therefore, the logotropic model is able to account 
both for cosmological and galactic observations remark¬ 
ably well. This may be a hint that DM and DE are the 
manifestation of a unique dark fluid. 

Finally, we have sketched some possible justifications 
of the logotropic equation of state in relation to the exis¬ 
tence of extra-dimensions (Cardassian model), in relation 
to Tsallis generalized thermodynamics, and in relation to 
scalar field theory and BECs. 

The fact that the Planck density pp enters in the lo¬ 
gotropic equation of state (|6| designed to model DM and 
DE is intriguing. It suggests that quantum mechanics 
manifests itself at the cosmological scale in relation to 
DE. This may be a hint for a fundamental theory of 
quantum gravity. This also suggests that the logotropic 
equation of state may be the limit of a more general equa¬ 
tion of state providing a possible unification of DE (p A ) 
in the late universe and inflation (vacuum energy pp) 
in the primordial universe. These open questions are a 


strong incentive to study the logotropic model further in 
future works. The phantom properties of the logotropic 
model will be discussed in a specific paper. 


Appendix A: Expression of the observational 
quantities in terms of the fundamental constants 


We enlight the remarkable feature that, in our the¬ 
ory, all the observational quantities can be predicted in 
terms of fundamental constants such as h, G, c, Hq, 
and Q mfi . We have B = l/[ln (87rc 5 /3f2 mi o?iGJ?o) — 1]- 
We introduce the notation % = [3f?(l — f2 m ,o)/2] = 

6.20 x 10" 2 . Then A = x 2 c 2 H'^/AttG ’= 3.06 x 
10 ~ 4 c 2 H 2 /G, E 0 = x^Hoc/AttG = 2.89 x 10~ 2 H 0 c/G, 
M h /r 2 h = X 0'(Zh)H o c/G = 4.30 x 10 ~ 2 H 0 c/G, 

v h/ r h = x9'(£h)H 0 c = 4.30 x 10~ 2 Hqc, v^/M b = 
Xd'(Zh)GH 0 c/f b = 4.30 x lO~ 2 GH 0 c/f b , a 0 = 

x0'(€h)H o c/f b = 4.30 x 10 ~ 2 H 0 c/f b , and M 300 /r 2 = 
xHqc/\/2G = 4.39 x 10~ 2 H 0 c/G. Noting that e A = 
Ac 2 /87 rG, e A = (1 - 0 )e 0 and Hq = (87rG/3c 2 )e 0 , 

we obtain A = 3(1 — Q m fi)H 2 = 2x 2 H 3 /B = 1.13 x 
10~ 35 s -2 . Therefore, the observational quantities can be 
expressed equivalently in terms of h, G, c, A and 


[1] F. Zwicky, Helv. Phys. Acta 6, 110 (1933) 

[2] V.C. Rubin, W.K. Ford, N. Thonnard, Astrophys. J. 
238, 471 (1980); M. Persic, P. Salucci, F. Stel, Mon. 
Not. R. astr. Soc. 281, 27 (1996) 

[3] A.G. Riess et al., Astron. J. 116, 1009 (1998); S. Perl- 
mutter et al., ApJ 517, 565 (1999); P. de Bernardis et 
al., Nature 404, 995 (2000); S. Hanany et al., ApJ 545, 
L5 (2000) 

[4] A. Einstein, Sitz. Konig. Preu. Akad. Wiss. 1, 142 (1917) 

[5] J.F. Navarro, C.S. Frenk, S.D.M. White, Mon. Not. R. 
astr. Soc. 462, 563 (1996) 

[6] A. Burkert, Astrophys. J. 447, L25 (1995) 

[7] G. Kauffmann, S.D.M. White, B. Guiderdoni, Mon. 
Not. R. astr. Soc. 264, 201 (1993); A. Klypin, A.V. 
Kravtsov, O. Valenzuela, Astrophys. J. 522, 82 (1999); 
M. Kamionkowski, A.R. Liddle, Phys. Rev. Lett. 84, 
4525 (2000) 

[8] F. Donato et al, Mon. Not. R. Astron. Soc. 397, 1169 
(2009) 

[9] R.B. Tully, J.R. Fisher, Astron. Astrophys. 54, 661 
(1977) 

[10] S.S. McGaugh, Astron. J. 143, 40 (2012) 

[11] L.E. Strigari et al, Nature 454, 1096 (2008) 

[12] H.J. de Vega, P. Salucci, N.G. Sanchez, Mon. Not. R. 
Astron. Soc. 442, 2717 (2014) 

[13] P.H. Chavanis, M. Lemou, F. Mehats, Phys. Rev. D 91, 
063531 (2015); P.H. Chavanis, M. Lemou, F. Mehats, 
arXiv: 1409.7840 

[14] A. Suarez, V.H. Robles, T. Matos, Astrophys. Space Sci. 
Proc. 38, 107 (2014); T. Rindler-Daller, P.R. Shapiro, 
Astrophys. Space Sci. Proc. 38, 163 (2014); P.H. Chava¬ 


nis, Self-gravitating Bose-Einstein condensates, in Quan¬ 
tum Aspects of Black Holes, edited by X. Calmet 
(Springer, 2015) 

[15] P. Steinhardt, in Critical Problems in Physics, edited by 
V.L. Fitch and D.R. Marlow (Princeton University Press, 
Princeton, NJ, 1997) 

[16] S. Weinberg, Rev. Mod. Phys. 61, 1 (1989); T. Padman- 
abhan, Phys. Rep. 380, 235 (2003) 

[17] G. Lemaitre, Proc. Nat. Acad. Sci. 20, 12 (1934); A.D. 
Sakharov, Dokl. Akad. Nauk SSSR 177, 70 (1967); Ya. 
B. Zeldovich, Sov. Phys. Uspek. 11, 381 (1968) 

[18] B. Ratra, J. Peebles, Phys. Rev. D 37, 321 (1988) 

[19] A. Kamenshchik, U. Moschella, V. Pasquier, Phys. Lett. 
B 511, 265 (2001) 

[20] H.B. Sandvik, M. Tegmark, M. Zaldarriaga, I. Waga, 
Phys. Rev. D 69, 123524 (2004); Z.H. Zhu, Astron. As¬ 
trophys. 423, 421 (2004) 

[21] P.H. Chavanis, C. Sire, Physica A 375, 140 (2007) 

[22] S. Weinberg, Gravitation and Cosmology (John Wiley, 

2002) 

[23] P.H. Chavanis, Eur. Phys. J. Plus 129, 222 (2014) 

[24] R.R. Caldwell, Phys. Lett. B 545, 23 (2002) 

[25] C. Destri, H.J. de Vega, N.G. Sanchez, New Astronomy 
22, 39 (2013) 

[26] S. Chandrasekhar, An Introduction to the Study of Stel¬ 
lar Structure (Dover, 1958) 

[27] A. Burkert, arXiv:1501.06604 (2015) 

[28] M. Milgrom, Astrophys. J. 270, 365 (1983) 

[29] K. Freese, M. Lewis, Phys. Lett. B 540, 1 (2002) 

[30] D. McLaughlin, R. Pudritz, Astrophys. J. 469, 194 
(1996) 







10 


[31] C. Tsallis, Introduction to Nonextensive Statistical Me¬ 
chanics (Springer, 2009) 

[32] P.H. Chavanis, Phys. Rev. D 84, 043531 (2011) 

[33] F. Dalfovo, S. Giorgini, L.P. Pitaevskii, S. Stringari, Rev. 
Mod. Phys. 71, 463 (1999) 


[34] E. Madelung, Zeit. F. Phys. 40, 322 (1927) 

[35] A. Suarez, P.H. Chavanis, Phys. Rev. D 92, 023510 
(2015) 



